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First measurement of the p spectral function in nuclear collisions 
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Abstract. The NA60 experiment at the CERN SPS has studied low-mass muon pairs in 158 AGeV In- 
In collisions. A strong excess of pairs is observed above the yield expected from neutral meson decays. 
The unprecedented sample size close to 400K events and the good mass resolution of about 2% made it 
possible to isolate the excess by subtraction of the decay sources. The shape of the resulting mass spectrum 
shows some non-trivial centrality dependence, but is largely consistent with a dominant contribution from 
tt^tt" ^ p ^ f^^ annihilation. The associated p spectral function exhibits considerable broadening, 
but essentially no shift in mass. The pr- differential mass spectra show the excess to be much stronger at 
low pt than at high pr- The results are compared to theoretical model predictions; they tend to rule out 
models linking hadron masses directly to the chiral condensate. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 12.38.Mh Quark-gluon plasma - 13.85.Qk Lepton Pairs 



1 Introduction 

Thermal dilepton production in the mass region < 1 GeV/c^ 
is largely mediated by the light vector mesons p, lu and 
(f>. Among these, the p(770 MeV/c^) is the most impor- 
tant, due to its strong coupling to the tttt channel and 
its short lifetime of only 1.3 fm/c, much shorter than the 
lifetime of the fireball. These properties have given it a 
key role as the test particle for "in-medium modifications" 
of hadron properties close to the QCD phase boundary. 
Changes both in width and in mass were originally sug- 
gested as precursor signatures of the chiral transition p. 
There seems to be some consensus now that the width of 
the p should increase towards the transition region, based 
on a number of quite different theoretical approaches ^ 
I21I3M1IS1- Oil the other hand, no consensus exists on how 
the mass of the p should change in approaching the transi- 
tion: predictions exist for a decrease PEEEIj a constant 
behavior and even an increase [S]. 

Experimentally, low-mass electron pair production was 
previously investigated at the CERN SPS by the CERES 
/NA45 experiment for p-Be/Au, S-Au and Pb-Au colli- 
sions [nillOllllj . The common feature of all results from 
nuclear collisions was an excess of the observed dilepton 
yield above the expected electromagnetic decays of neutral 
mesons, by a factor of 2-3, for masses above 0.2 GeV/c^. 
The surplus yield has generally been attributed to direct 
thermal radiation from the fireball, dominated by pion 
annihilation tt~^tt^ —tp—t l^l^ with an intermediate p 
which is strongly modified by the medium. Statistical ac- 
curacy and mass resolution of the data were, however, 
not sufficient to reach the sensitivity required to assess in 
detail the character of the in-medium changes. The new 



experiment NA60 at the CERN SPS has now achieved a 
decisive breakthrough in this field. 



2 Apparatus and data analysis 

The apparatus is based on the muon spectrometer pre- 
viously used by NA50, and a newly added telescope of 
radiation-tolerant silicon pixel detectors, embedded inside 
a 2.5 T dipolc magnet in the vertex region (12t.l3; . Match- 
ing of the muon tracks before and after the hadron ab- 
sorber, both in angular and momentum space, improves 
the dimuon mass resolution in the region of the light vec- 
tor mesons from ~80 to ~20 MeV/c^ and also decreases 
the combinatorial background of muons from tt and K de- 
cays. Moreover, the additional bend by the dipole field 
leads to a strong increase of the detector acceptance for 
opposite-sign dimuons of low mass and low transverse mo- 
mentum. The rapidity coverage is 3.3<y<4.3 for the p, at 
low Pt (compared to 3<y<4 for the J/V')- Finally, the se- 
lective dimuon trigger and the fast readout speed of the 
pixel telescope allow the experiment to run at very high 
luminosities. 

The results reported here were obtained from the anal- 
ysis of data taken in 2003 with a 158 AGcV indium beam, 
incident on a segmented indium target of seven disks with 
a total of 18% (In-In) interaction length. At an average 
beam intensity of 5T0^ ions per 5 s burst, about 3-10^^ 
ions were delivered to the experiment, and a total of 230 
million dimuon triggers were recorded on tape. The data 
reconstruction starts with the muon-spectrometer tracks. 
Next, pattern recognition and tracking in the vertex tele- 
scope are done; the interaction vertex in the target is 
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reconstructed with a resolution of ^200 fim for the z- 
coordinate and 10-20 /im in the transverse plane. Only 
events with one vertex are kept; interaction pileup and 
reinteractions of secondaries and fragments are thus re- 
jected. Finally, each muon-spectrometer track is extrapo- 
lated to the vertex region and matched to the tracks from 
the vertex telescope. 

The combinatorial background of uncorrelated muon 
pairs mainly originating from tt and K decays is deter- 
mined using a mixed-event technique '14'. Two single mu- 
ons from different like-sign dimuon triggers are combined 
into muon pairs in such a way as to accurately account for 
details of the acceptance and trigger conditions. The qual- 
ity of the mixed-event technique can be judged by com- 
paring the like-sign distributions generated from mixed 
events with the measured like-sign distributions. It is re- 
markable that the two agree to within ^1% over a dy- 
namic range of 4 orders of magnitude in the steeply falling 
mass spectrum jl4| . After subtraction of the combinato- 
rial background, the remaining opposite-sign pairs still 
contain "signal" fake matches, i.e. associations of muons 
to non-muon tracks in the pixel telescope. This contribu- 
tion is only 7% of the combinatorial background level. It 
has been determined in the present analysis by an over- 
lay Monte Carlo method. We have verified that an event- 
mixing technique gives the same results, both in shape 
and in yield, within better than 5%. More details on the 
experimental apparatus and data analysis will be given in 
a forthcoming extended paper; for now see |14lll5j . 



3 Results 

A significant part of the results presented in this paper 
has recently been published ^B]- That part will therefore 
be less extensively treated than other parts published here 
for the first time. 

Fig. n shows the opposite-sign, background and sig- 
nal dimuon mass spectra, integrated over all collision cen- 
tralities. After subtracting the combinatorial background 
and the signal fake matches, the resulting net spectrum 
contains about 360 000 muon pairs in the mass range 0- 
2 GeV/c^ of Fig. H] roughly 50% of the total available 
statistics. The average charged-particle multiplicity den- 
sity measured by the vertex tracker is dNch/dr] =120, the 
average signal-to-background ratio is 1/7. For the first 
time in nuclear collisions, the vector mesons lo and (j) 
are completely resolved in the dilepton channel; even the 
r]-^fj,fj, decay is seen. The mass resolution at the uj is 20 
MeV/c^. The subsequent analysis is done in four classes 
of collision centrality defined through the charged-particle 
multiplicity density: peripheral (4-30), semiperipheral (30- 
110), semicentral (110-170) and central (170-240). The 
signal-to-background ratios associated with the individual 
classes are 2, 1/3, 1/8 and 1/11, respectively. 

The peripheral data can essentially be described by the 
expected electromagnetic decays of the neutral mesons. 
Muon pairs produced from the 2-body decays of the j], p, 
LO and (j) resonances and the Dalitz decays of the ry, ry and 
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Fig. 1. Mass spectra of the opposite-sign dimuons (upper 
histogram), combinatorial background (dashed), signal fake 
matches (dashed-dotted), and resulting signal (histogram with 
error bars). 



Lu were simulated using the improved hadron decay gen- 
erator GENESIS [HI, while GEANT was used for trans- 
port through the detectors. Four free parameters (apart 
from the overall normalization) were used in the fit of this 
"hadron decay cocktail" to the peripheral data: the cross 
section ratios jy/w, p/to and '/>/a;, and the level of D meson 
pair decays; the ratio rj /rj was kept fixed at 0.12 |nill7j . 
The fits were done without pt selection, but also indepen- 
dently in three windows of dimuon transverse momentum: 
Pt<0.5, 0.5<pt<1 and pt>1 GeV/c. Data and fits, in- 
cluding an illustration of the individual sources, are shown 
in Fig. 12 for all pr (upper) and for the particular selection 
Pt<0.5 GeV/c (lower). The fit quahty is good throughout, 
even in the critical acceptance-suppressed 77-Dalitz region 
at low mass and low p-p. 

The quantitative fit results in terms of the cross section 
ratios, corrected for acceptance and extrapolated to full 
phase space (meaning here full ranges in y and pr) are dis- 
played in Fig. 01 (upper), where the horizontal lines (with 
the label <>) indicate the fit values without p^ selection, 
and the data points (with statistical errors) the values 
obtained from the three different pr windows. The sys- 
tematic errors are of order 10% in all cases, dominated by 
those of the branching ratios. The rj/uj ratio agrees, within 
<10%, with the literature average for p-p, p-Be 0. The 
(f>/ijj ratio is higher than the p-p, p-Be average, reflecting 
some (j) enhancement already in peripheral nuclear colli- 
sions. Both ratios arc, within 10%, independent of the pair 
pr. This implies that the GENESIS input assumptions 
(used in the extrapolation to full pt) as well as the accep- 
tance corrections vs. pt are correct on the level of 10%. 
As shown in Fig. (lower), the acceptance variations with 
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Fig. 2. Fits of hadron decay cocktail to the peripheral data 
for all pt (upper) and pt<0.5 GeV/c (lower), showing also the 
individual contributions. 
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Fig. 3. Upper: Particle cross-section ratios for all pt (<>) and 
for three pr windows, extrapolated to full phase space. Lower: 
NA60 acceptance relative to 47r for different mass windows. 



Pt are minor for the to and the (f), but very strong, over 
two orders of magnitude, for the ry Dahtz mode (M<0.4 
GeV/c^). The accuracy level of 10% reached in under- 
standing the acceptance is therefore truly remarkable. 

The particle ratio p/uj behaves in a different way rel- 
ative to the other two. It decreases with p^, but remains 
significantly higher than the p-p, p-Be average ;9 through- 
out. This suggests that some tttt annihilation, enhancing 
the yield of the p, contributes already in peripheral colli- 
sions (see below). 

In the more central bins, a fit procedure is ruled out, 
due to the existence of a strong excess with a priori un- 
known characteristics. We have therefore used a novel pro- 
cedure as shown in Fig. 01 made possible by the high data 
quality. The excess is isolated by subtracting the cocktail, 
without the p, from the data. The cocktail is fixed, sep- 
arately for the major sources and in each centrality bin, 
by a "conservative" approach. The yields of the narrow 
vector mesons co and are fixed so as to get, after sub- 
traction, a smooth underlying continuum. For the ry, an 
upper limit is defined by "saturating" the measured data 
in the region close to 0.2 GeV/c^; this implies the excess 



to vanish at very low mass, by construction. The r/ reso- 
nance and uj Dalitz decays are now bound as well; r/ /rj is 
fixed as before. The cocktail p (shown in Figs.|SlEllHland|51 
for illustration purposes) is bound by the ratio p/cu— 1.2. 
The accuracy in the determination of the to and yields 
by this subtraction procedure is on the level of very few 
%, due to the remarkable local sensitivity, and not much 
worse for the 77. 

The excess mass spectra for all 4 multiplicity bins, re- 
sulting from subtraction of the "conservative" hadron de- 
cay cocktail from the measured data, are shown in Figs. El 
andOfor all pr and the particular selection pr<0.5 GeV/c, 
respectively. The cocktail p and the level of charm decays, 
found in the three upper centrality bins to be about 1/3 
of the measured yield in the mass interval 1.2<M<1.4 
GeV/c^ are shown for comparison. The qualitative 
features of the spectra are striking: a peaked structure is 
seen in all cases, broadening strongly with centrality, but 
remaining essentially centered around the position of the 
nominal p pole. At the same time, the total yield increases 
relative to the cocktail p, their ratio (for M<0.9 GeV/c^) 
reaching values of 4 for all pt, even close to 8 for pT<0.5 
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Fig. 4. Isolation of an excess above the hadron decay cocktail 
(see text). Total data (open circles), individual cocktail sources 
(solid), difference data (thick triangles), sum of cocktail sources 
and difference data (dashed). 



GeV/c, in the most central bin. Such values are consis- 
tent with the CERES ^T] results, if the latter are also 
referred to the cocktail p and rescaled to In-In. The er- 
rors shown are purely statistical. The systematic errors are 
dominantly connected to the uncertainties in the level of 
the combinatorial background, less so to the fake matches. 
For the data without pT selection, they are estimated to 
be about 3%, 12%, 25% and 25% for the 4 centralities 
in the broad continuum region, while the /o-like structure 
above the continuum is much more robust. Uncertainties 
associated with the subtraction of the hadron decay cock- 
tail reach locally up to 15%, dominated by those of the 
w-Dalitz form factor. For the data with pT<0.5 GeV/c, 
the systematic errors are still under investigation. 



For the data without p^ selection, a quantitative anal- 
ysis of the shape of the excess mass spectra vs. centrality 
has been performed, using a finer subdivision into 12 cen- 
trality bins. Referring to Fig. El the data were subdivided 
into three mass windows with equal widths: 0.44<M<0.64 
(L^Lower), 0.64<M<0.84 (C=Center), and 0.84<M<1.04 
GeV/c^ (U=Upper). From the yields in these windows, 
a peak yield R = C - 1/2(L-|-U) and a continuum yield 
3/2 (L+U) can be defined. Fig. |7| (upper) shows the ra- 
tios of peak/p (R/p), continuum/p (3/2(L-|-U)/p), and 
peak/continuum (RR), where p stands for the cocktail p. 
The errors shown are purely statistical. The relative sys- 
tematic errors between neighboring points are small com- 
pared to the statistical errors and can therefore be ignored, 
while the common systematic errors are related to those 
discussed in the previous paragraph and are of no rele- 
vance here. The ratio peak/p is seen to decrease from the 
most peripheral to the most central bin by nearly a factor 
of 2, ruling out the naive view that the shape can simply 
be explained by the cocktail p residing on a broad con- 
tinuum, independent of centrality. The ratio continuum/p 
shows a fast initial rise, followed by a more flat and then 
another more rapid rise beyond dNc/i/dy = 100; this be- 
havior is statistically significant. The sum of the two ratios 
is the total enhancement factor relative to the cocktail p; 
it reaches about 5.5 in the most central bin. The ratio of 
the two ratios, peak/continuum, amplifies the two sepa- 
rate tendencies: a fast decay, a nearly constant part, and 
a decline by a further factor of 2 beyond dNc/i/dy =100, 
though with large errors due to R. 

A completely independent shape analysis was done by 
just evaluating the RMS = a/{M2) - (M)^ of the mass 
spectra in the single total mass interval 0.44<M<1.04 
GeV/c^. The results are shown in Fig. |7| (lower), both 
for the full data and after subtraction of charm on the 
level discussed for Figs. El El (except for the most periph- 
eral bin where the continuum yield above 1 GeV/c^ was 
assumed to be 100% charm). They are in perfect qualita- 
tive agreement with the more "microscopic" shape analy- 
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Fig. 5. Excess mass spectra of dimuons without pt selection. 
The cocktail p (sohd) and the level of uncorrelated charm de- 
cays (dashed) are shown for comparison. For errors see text. 
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Fig. 6. Excess mass spectra of dimuons for pt<0.5 GeV/c. 
The cocktail p (solid) and the level of uncorrelated charm de- 
cays (dashed) are shown for comparison. For errors see text. 
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sis discussed before, rising from values close to the cocktail 
p all the way up to nearly a flat-continuum value. The ex- 
tra rise beyond dNc/i/dy — 100 is highly significant here, 
due to the very small statistical and systematic errors. 



4 Comparison to theoretical models 

The qualitative features of the excess mass spectra shown 
in Figs. 13 and are consistent with an interpretation as 
direct thermal radiation from the fireball, dominated by 
TTTT annihilation. A quantitative comparison of the data to 
the respective theoretical models can either be done at the 
input of the experiment, requiring acceptance correction 
of the data, or at the output, requiring propagation of the 
theoretical results through the experimental acceptance. 
All the data contained in this paper have so far not been 
corrected for acceptance, and therefore only the second al- 
ternative is available at present (the first one being under 
preparation). To help intuition, Fig. |H1 iflustrates the ef- 
fects of acceptance propagation for the particularly trans- 
parent case of qq annihilation, associated with a uniform 



spectral function jl8| . By coincidence, without p^ selec- 
tion, the resulting mass spectrum at the output is also 
uniform within 10% up to about 1.0 GeV/c^, resembling 
the shape of the spectral function at the input. In other 
words, the always existing steep rise of the theoretical in- 
put at low masses (Fig. ISJ , due to the photon propagator 
and a Boltzmann-like factor 6 ,7 , is just about com- 
pensated by the falling acceptance in this region as long as 
no pt cut is applied. Variations of the input p^ spectrum 
within reasonable physics limits affect the fiatness of the 
output by at most 20%. Strong cuts in p^ like <0.5 or >1 
GeV/c, however, completely invalidate the argument. 

On the basis of this discussion, the excess mass spectra 
shown in Fig. can approximately be interpreted as spec- 
tral functions of the p, averaged over momenta and the 
complete space-time evolution of the fireball. The broad 
continuum- like part of the spectra may then reflect the 
early history close to the QCD boundary with a nearly 
divergent width, while the narrow peak on top may just 
be due to the late part close to thermal freeze-out, ap- 
proaching the nominal width. The p-r-cut data shown in 
Fig. 1^1 on the other hand, do not allow such interpreta- 
tion, due to the extreme acceptance cut on the low-mass 
side of the p; compare also Fig. 13 (lower). 

Among the many different theoretical predictions for 
the properties of the intermediate p mentioned in the in- 
troduction, only few have been brought to a level suitable 
for a quantitative comparison to the data. Before the first 
release of the present data in 2005, the in-medium broad- 
ening scenario of and the moving mass scenario re- 
lated to inilZI were evaluated for In-In at dNch/dr7=140, 
using exactly the same fireball evolution model extrapo- 
lated from Pb-Pb, and taking explicit account of temper- 
ature as well as of baryon density ^H]. In Fig. ^ these 
predictions (as well as the unmodified p) are confronted 
with the data for the semicentral bin. The integrals of 
the theoretical spectra are independently normalized to 
the data in the mass interval M<0.9 GeV/c^ in order to 
concentrate on the spectral shapes, independent of the un- 
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Fig. 8. Propagation of thermal qq radiation, based on a uni- 
form spectral function, through the NA60 acceptance without 
any pr selection. The resulting spectrum is, by coincidence, 
also uniform. 
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(which should also be viewed as a systematic error). 



certainties of the fireball evolution. Note again that data 
and predictions can only be interpreted as space-time av- 
eraged spectral functions for the part without pr selection 
(Fig. 121 center), but not for the other two. Irrespective of 
the choice of the pr window, however, some general con- 
clusions can be drawn. The unmodified p is clearly ruled 
out. The moving mass scenario related to Brown/Rho scal- 
ing, which fit the CERES data |KIITT] . is also ruled out, 
showing the much improved discrimination power of the 
present data. Only the broadening scenario appears to be 
in fair agreement with the data. 

In the meantime, the release of the data in 2005 has 
triggered a number of new theoretical developments. Con- 
trary to initial critics 19 , Brown/Rho scaling could not 
be saved by varying fireball and other parameters within 
extremes, including switching out the effects of tempera- 
ture altogether 20,21 . The excess of the data at M>0.9 
GeV/c^ may be related to the prompt dimuon excess found 
by NA60 in the intermediate mass region It is not 
accounted for by the model results shown in Fig. and 
it is presently (nearly quantitatively) described either by 
hadronic processes like 47r, Btt.. (including vector-axial- 
vector mixing) |22|. or by partonic processes like qq anni- 
hilation 1^ . This is a challenging theoretical ambiguity to 
be solved in the future. A chiral virial approach has also 
been able to nearly quantitatively describe the data (24j. 



5 Conclusions 

The data unambiguously show that the p primarily broad- 
ens in In-In collisions, but does not show any shift in mass. 
Consequently, model comparisons favor broadening sce- 
narios, but tend to rule out moving-mass scenarios cou- 
pled directly to the chiral condensate. The issue of vector- 
axialvector mixing, also sensitive to chiral restoration, re- 
mains somewhat open at present. We expect that precise 
Pt dependences, presently under investigation, will give 
more insight into the different sources operating in differ- 
ent mass regions. 
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